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Abstract. We present the kinematics of the molecular gas in the barred spiral galaxy M83 (NGC5236). The 
study is based on ^^CO(J=1-0 and 2-1) observations with the Swedish-ESO Submillimetre Telescope (SEST). 
Iso- velocity maps of the entire optical disk, 10' x 10' or 13 x 13kpc, are produced. They show the pattern of an 
inclined, rotating disk, but also the effects of streaming motions along the spiral arms. A dynamical mass of 
' about 6 X 10^" Mq is estimated by fitting the rotation curve of an exponential disk model to these data. The 

I gas constitutes about 13% of the disk mass. The pattern speed is determined from the residual velocity pattern 

, The locations of various resonances are discussed. The molecular gas velocity dispersion is determined, and a 

' trend of decreasing dispersion with increasing galactocentric radius is found. A total gas (H[2+H[l+He) mass 

I surface density map is presented, and compared to the critical density for star formation of an isothermal gaseous 

. disk. The critical density is exceeded in the spiral arms, but not in the interarm regions. The locations of Giant 

' Molecular Associations (GMAs) and Hll regions are consistent with this scenario of dynamically induced star 

formation. 
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1. Introduction lines in the optical regime that can be used to trace the 

motion of the disk. However, since the kinematics of these 

The kinematics of the gas in a galaxy provides fundamen- regions are affected by the starformation process it can be 

5^ ; tal clues to its dynamics and evolution. An excellent probe difficult to interpret such data. Also, the high extinction 
of the kinematics of disk galaxies is the molecular gas (as starforming regions further complicates matter, 
traced by e.g. CO radio line emission). Since the molecu- ... 
, • IT • 1 ii J J.- J J J The molecular gas kinematics is also an important 
lar gas is coUisional, the random motions are reduced and i n i r a i i 

.1, 1 J ID fii,- 4.U 1 probe of the star formation process. A spiral density wave 

the disk IS dynamically cooled, l^ecause of this the molec- , , . r . ■ , 

, • ui i J T ii -11 J i 1 i IS an efficient mechanism for concentrating atomic and 

ular gas is able to respond, both quickly and strongly, to ° 

, • J.U -4- 4.- 1 i 4.- 1 Ai ■ u J molecular hydrogen m spiral features, but m order to ini- 

cfianges m the gravitational potential. Atomic fiydrogen r 

, , 1,- 1, • fi J u 4- ■ 4.V. 1 tiate the star formation, the disk also needs to be com- 

is also a tracer which is often used, but since the scale ' 

1 . 1 , f- XT, • 1 4-- 1 1 J 4- 4-1, 4- f 4- V, pressed m the vertical direction. In this smaller-scale pro- 
height of Hi is relatively large compared to that oi the ^ ^ 

, , -4. • 1 -i- 4. 4.1, J ■ ■ 4-u cess, dissipative processes like cloud-cloud collisions, m- 

moiecular gas it is less sensitive to the dynamics m the ii i nii ii 

,.1 , a-4-jur 4- CJ4.4.UJ-1 teraction with the galactic magnetic field, and the Jeans 

disk, and more afiected by forces not confined to the disk . , .,. i i f ■ i i • 

, / 1 , r 1 • \ A 1 instability help enhancing the density. Once a certain 

plane (e.g., shock waves from supernova explosions). Also, , , , . . , , . r ■ ■ 

i r 1 • 11 1 1-4.4-1 TT • threshold density is reached, massive star formation is 

the centers of galaxies usually show very little H i emis- , , , 

Q, n • • ■ ■ 4- u -(• ■ ■ initiated since the velocity dispersion no l onger can neu- 

sion. btarformmg regions give rise to a number of emission '^^ i <^ ii-^ — i 

tralize the self-g ravity of the gaseous disk l|Toomrdll964i 



traiize tne selt-g r 
lKennicuttlll989|i . 



Send offprint requests to: Andreas Andersson Lundgren (an- 

dreas@astro.su.se) Two large surveys of the kinematics of spiral galax- 

* Based on observations collected with the Swedish- ies using CO radio line emission have recently been com- 

ESO Submillimetre Telescope at the European Southern plet ed: "B IMA SONG" done with the BIMA interferome- 



Observatory, La Silla, Chile ter l)Regan et al...200L: .Heifer et al...200,'^) and "The Virgo 
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HiKh-Resolut ion CO Survey" do ne with the Nobeyama in- 
terferometer ijSofue et alJl2003j) . Both surveys noted that 
some galaxies have a "twin-peaks" distribution of the 
molecular gas in the center, and it is suggested that this 
is a result of molecular gas agglomeration at the inner 
Lindblad resonance. More results on the kinematics will 
come in subsequent papers from these groups. Many other 
galaxies have been observed in CO with interferometers 
and among the most detailed studies of the global kine- 
matics are studies done on M 100 w ith the IRAM interfer- 
ometer ijCarcia-Burillo et al.lll998|) and one on M51 with 
the OVRO interferometer ijAalto et aLlll999() . In M 100 
the kinematics of the molecular gas was modeled using a 
hydrodynamic code, and it was shown that the stellar bar 
and the spiral pattern have a common pattern speed, while 
the nuclear bar is rotating significantly faster. In M 51 the 
authors found streaming motions of very high velocities 
(60-150 km s~^) and that two pattern speeds coexist - 
one main mode and a second due to the interaction with 
NGC5195. 

Detailed studies done with single-dish telescopes of 
the dynamics of nearby galaxies with low inclinations 
are rare. The most well-studied galaxy is M51, which 
has been observ ed in CO radio lin e emission with 
the Onsala 20m llRvdbeck et aLlHiil . the IRAM 30m 
( Garcia-B urillo et'^dll99 3^, and the Nobeyarn a 45m tele- 
scope l|Kuno et al.lll995t iKuno fc Nakailll997(l . The last 
group showed that the gaseous disk is gravitationally un- 
stable in the arms, but not in the inter-arm regions, and 
that the deviations from circular motion seen in the gas 
are in accordance with density-wave theory. 

The kinematics of M 83 has been observed using many 
tracers. In particular, CO radio line d ata results have 
been presented by several authors; see iLunderen et alJ 
I 2004ll rPaper H for a sum mary of CO observations of 
M 83. ICombes et all l(l978l) produced the first CO ro- 



tation curve based on data in 7 positions. Position- 
velocity diagrams of the nuclear region a long the mi- 
nor an d major axes have bee n publis hed bv [ Ha nda et all 
lll990t). Ilsrael fc BaasI (l200ll) . and .Dumke et ah (.2001,). 
IPetitpas fc WilsonI l|l998() "presented channel maps of the 
central 40" x 40", and discus sed inner Lindbla. d resonances 
and gas flow along the bar. iLord fc KennevI (^91j) used 
the OVRO interferometer to observe the eastern arm of 
M83 in the CO(J=1~0) line. They presented channel 
maps and discussed the relation between CO, Ha, and 
H I. The same authors also observed the western bar end 
i Kennev fc Lordlll99ll) . and they d iscussed, among other 
things, orbit crowding. iRand et"!.. (,1999.1 revisited the 
eastern arm, extended the map, and discussed stream- 
ing motions and disk stability. The most extensive work 
on the molecular gas kinematics prior to this paper was 
done by Crosthwaite et al. (2002), who mapped the disk 
in CO(J=1-0 and J=2-l) and presented a CO rotation 
curve, iso-velocity patterns, and velocity dispersion pat- 
terns covering the optical disk. 

The velocity pattern in Hi, both in the disk 
and in the Hi envelope, and discussions of an 



Hi warp have been prese nted in paper s base d on 
both interferomet er data llRoestad et alJ Il974l) an d 
single-dish data fanchtmeier fc Bohnenstengell Il98ll) . 
iTilanus fc AllenI ( 19931) compared velocitv information ob- 
tained from H I and H/3 data, and discussed shocks, large 
scale H2 dissociation, and the relation between the dust 
an d Hi. 

IComtel lITflsl and Ide Va,ucon1eurs et a,1.l llT9si pre- 
sented Fabry-Perot Ha observations of M83. The latter 
group also discussed the pattern speed, the mass distribu- 
tion, and the residual velocities. 

In this paper, we present the kinematics of the molec- 
ular gas, as seen in the C0( J=l-0) and C0( J=2-l) lines, 
in the galaxy M 83. Two-dimensional velocity distribution 
and velocity dispersion maps are shown, and a rotation 
curve based on the assumption of an exponential disk mass 
distribution is derived. By subtracting the circular veloci- 
ties from the observed data we get the residual velocities. 
Lastly, the location and formation of GMAs and the sta- 
bility of the gaseous disk are discussed. 

2. Observations 

Our CO(J=1-0 and 2-1) observations were presented 
in Paper I. The data cover the entire optical disk, and 
is sampled at intervals of 7"-ll". The spacing is small 
enough to facilitate the use of a MEM-deconvolution 
technique. The data grid is oriented along the equato- 
rial coordinate system and centered on the coordinates 
13' |36°'59^4, -29°52'05" which is the optical center given 
bv lde Vaucouleurs et al.l l(l97fi|) . Since the grid spacing is 
small relative to the beam size, we can gain in S/N by 
convolving the raw data sets, without significantly affect- 
ing the angular resolution. The results in this paper are 
therefore based on convolved CO(J=1-0) data (convolu- 
tion kernel FWHM 20", final angular resolution 49"), con- 
volved C0(J=2-1) data (kernel 15", resolution 27") and 
two MEM-deconvolved CO(J=1-0) and C0(J=2-1) data 
sets with angular resolutions of 22" and 14", respectively. 
In order to compare data at the same angular resolution, 
data cubes with angular resolutions of 49" (based on the 
raw data sets) and 27" (based on the MEM-deconvolved 
data sets) were produced. The C0( J=2-l) data have been 
corrected for emission in the error beam. For a detailed 
description of the data, see Paper I. 

3. Kinematics 

3.1. Channel maps 

In Fig. n channel maps ranging from 400 to 620 km s~^ of 
the MEM-deconvolved C0( J=l-0) data cube, where each 
map covers a velocity range of 15 km s^^, are shown. They 
show the expected overall pattern of an inclined, rotating 
disk. In addition, the arms are clearly resolved. On the 
global scale the most notable deviation from the expected 
is the outer arm on the NW side, which seems slightly 
offset from the expected velocity, an indication of stream- 
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Table 1. General Parameters of M83. 



Morphological Type" 
IR center:'' 
R.A. (J2000) 
Decl. (J2000) 

LSR systemic velocity (opt)'' 
Distance'* 
Position angle'' 
Inclination-^ 

Holmberg diameter (Do)" 

Mhi" 



SAB(s)c 

13''37™00!8 

-29''51'56" 

506 km s"^ 

4.5 Mpc 

45° 

24° 

14^6 

7.7x10'* Mo 
3.9x10^ Mo 



" ^^V^KOuleurs et al^ ^^7^; 

(1994); ° Comtc 71981) : 
Huchtmcior & Bohncnstengel 
'Lundgren et al. (2004) 



^ofu^^^W^camateu 
Thim et al. (2003); " 



T 



Talbot et al. (1979 



ing motions or a warp. The prominent central region is 
visible throughout the entire velocity range. The emission 
peak moves from (+9",+2") to (-4",-8") relative to the 
IR center (See Tabled]) - a shift of 21", which corresponds 
to 460 pc. These positions match with in a few arcseconds 
the tw o nuclear components found bv lPetitoas fc WilsonI 
l|l998tl in their CO(J=3-2 and 4-3) data. 

In Fig. 121 the same channel maps in the MEM- 
deconvolved C0(J=2-1) data cube are shown. In order to 
highlight the finer details provided by the better resolu- 
tion in this data set, only the inner 140" x 140" are shown. 
Here, the two central components are clearly resolved. 

3.2. Rotation curve and kinematic fit 

There exist many proced ures for measur i ng th e velocity 
of a spectral feature (see ISofue fc RubinI l)200l|) for a re- 
cent review of methods). We chose the intensity- weighted- 
velocity method to be able to directly compare with results 
obtained using the same method applied to other velocity 
tracers, such as Ha and Hi. For a galaxy with a low incli- 
nation, like M 83, most methods give almost identical rota- 
tion curves. However, in the nuclear region, where the ve- 
locity changes very fast with respect to the beam size, it is 
likely that the intensity-weighted-velocity method under- 
estimates the true rotational velocity. In Section Hi 5. II we 
also derive the rotation curve using the envelope-tracing 
method ( Sofuc et al. 1997). 

In each spectrum the "average" velocity is estimated 
by calculating the intensity-weighted velocity. In order 
to minimize the effect of baseline variations, a sliding- 
window technique was used, where at each position in- 
tegration was performed inside a velocity window whose 
center and width depend on the position within the galaxy. 
We checked carefully that the window always covered all 
of the emission. 

To obtain the best-fit circular velocities in the galaxy, 
the kinematic data have to be compensated for inclination, 
position angle, systemic velocity, and kinematic-center off- 
set. In order to estimate the values of these parameters. 



a rotation curve is introduced to which the velocity data 
sets can be compared. Here we have chosen the theoretical 
rotation curve produced by an exponential disk potential. 
The circular velocities in such a disk are given by: 



y 



vi{y) = 2GMd^ [h{y)Ko{y) - h{y)Ki{y)] 



(1) 



where Md is the total mass of the exp onential disk, Rd the 
disk scale length, and y = R/2Rd l|Binnev & Tremain3 
Il987l) . /„ and i^n are the modified Bessel functions of the 
first and second kind, respectively. 

We expect that the gas will not move in circular orbits 
in the bar region. Therefore, data obtained within 1' of the 
nucleus were excluded in the fit. In order to minimize the 
influence of radial motions data from areas further away 
than 40'' from the major axis were also excluded. With 
this method it is not possible to determine the inclination 
angle, since the unique effects of variations of this parame- 
ter are seen as a changing opening angle of the iso- velocity 
contours along the minor axis^. We adopted an inclination 
angle of 24°, which is in agreement wit h both the H i iso- 
velocity contours and optical isophotes llComtelll98lh . All 
other parameters were fitted simultaneously to the kine- 
matic data sets obtained from the convolved CO(J=1-0) 
and C0( J=2-l) data cubes. The fitting error s were mini- 
mized using a Levenberg-Marquardt method ijPress et alJ 
[l99i). We found that the kinematic center lies about 12" 
SE of the IR center (see Table ^ , but within the errors 
of the fit the two positions agree. The rotation curves de- 
rived for the south-western and north-eastern quadrants 
separately differed markedly in the inner 1'. This discrep- 
ancy disappeared when the kinematic offset parameters 
were fixed to match the location of the IR center. We 
fitted the remaining parameters, and the best-fit values 
(with errors) are given in Table El The same procedure 
was applied to the MEM-deconvolved data sets, and the 
same results were obtained (within the errors). Our re- 
sul ts agree closely with th ose obtained from the Hi data 
bv lTilanus fc Allen! l)l993(l (hereafter TA). 

The rotati on curve fit results in a dynamical mass of 
6 X 10^° M0. ((Crosthwaite et al.l l|2002l) derived a dynam- 
ical mass of 8 X 10^'^ Mq (for a distance of 4.5 Mpc) using 
a Brandt model rotation curve.) The total mass of molec- 
ular and atomic hydrogen gas within the Holmberg limit 
(5.5 X 10^ M0, see Paper I) constitutes about 9%±2% 
of the estimated disk mass. Including helium the total 
gas mass is about 13% of the disk mass. This fraction de- 
pends somewhat on the choice of the inclination angle e.g., 
changing the inclination from 24 to 26 degrees increases 
the hydrogen gas fraction from 9% to 10%. 

The circular velocities rise linearly (and slightly slower 
than the model rotation curve) in the central regions and 
level off at larger radii, upper panel in Fig. O In this fig- 
ure we have f or c omparison also plotted circular veloci- 
ties from Ha i^Vaucoulcurs et al. 1983) and Hi (TA). 

^ Along the major axis, an increase of inclination results in a 
higher rotation curve, which can also be interpreted as a larger 
disk mass - thus a degenerate problem. 
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Fig. 1. Channel maps of the MEM-deconvolved CO(J=1-0) data set ranging from 400 to 640 km s~^ with a spacing 
of 15 km s^^. The velocity interval is indicated in the upper left corner of each panel. Contour levels are 0.5, 1.5, 3, 
6, 12, and 24 K km s~^. The cross marks the IR center. 



Inspection of position-velocity diagrams along the major 
axis (Fig.O indicates that there are substantial deviations 
from the rotation curve due to e.g. streaming motions (see 
Sect. El. 

Although not readily visible from the azimuthally- 
averaged rotation curve, a central component of relatively 
low mass may exist. Indeed, the position-velocity dia- 
grams along the major axis, shown in Fig. [SJ suggest 
a more rapidly rising rotation curve with a first maxi- 
mum at « 15". Adopting, as an example (and consistent 



with our estimate from the position-velocity diagrams in 
Sect. 13.5. l|l . an inner (nuclear) exponential disk^ with a 
mass of 3 X 10^ Mq, and a disk scale length of 50 pc, and 
plotting the new rotation curve together with the circu- 
lar velocities obtained from the deconvolved C0(J=2-1) 
data set (lower panel in Fig.lHJ, it is clear that the effects 



^ The type of mass distribution has relatively little effect on 
the behavior of the rotation curve as long as the scale length 
of the mass distribution is sufficiently small. 
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Fig. 2. Channel maps of the central parts of the MEM-deconvolved C0(J=2-1) data set. The velocity interval and 
contour levels are the same as in Fig. ^ The cross marks the IR center. 

Table 2. Best-fit parameters for our rotating-disk model. The kinematic center was kept fixed to the IR center 
(coordinates given in Table The data used are those of the convolved data sets. 





CO(J 


=1-0) 


CO(J 


-2-1) 




Value 


Error (la) 


Value 


Error (la) 


Kinematic center 


IR center 


fixed 


IR center 


fixed 


Position angle 


46° 


5° 


46° 


5° 


Inclination 


24° 


fixed 


24° 


fixed 


Systemic vel. (LSR) [km s"^] 


511.8 


0.6 


511.5 


0.7 


Disk mass [10^0 Mq] 


5.9 


0.7 


6.2 


0.8 


Disk scale length [kpc] 


2.7 


0.2 


2.9 


0.2 
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ftp (corotation, CR) or when 
nm(Cl — Qp) = ±K, where m is the number of arms (in 



of such a mass distribution can very well be hidden in our 
data. (The mass and the scale length of the outer disk 
remain the same.) 



3.3. Resonances 

Objects orbiting the center of a galaxy oscillates radially 
in addition to their motion in circular or nearly-circular 
orbits. If the frequency of these radial oscillations, called 
the epicycle frequency (k), is commensurate with the an- 
gular velocity of the bar and spiral arms, resonances oc- 
cur. Depending on the conditions, the resonance can trap 
objects in a density wave or deplete a region. These res- 
onances occur when 

n 

this case 2), n an integer number > 1, f2 the angular ve- 
locity of the gas and stars, and flp the pattern speed. The 
case n—1 results in the Lindblad resonances, while n = 2 
corresponds to the lowest order of the ultraharmonic reso- 
nances. In barred galaxies it is generally assumed that the 
bar and the spiral pattern rotate with a common pattern 
speed. In the following discussion, we will assume that the 
bar and the spiral arms have a common pattern speed of 
47 km s^i kpc~i (see Sect.lX^. 

This results in a corotation radius (CR) of «170" (3.7 
kpc), an outer UHR (oUHR) at «240" (5.2 kpc), and 
an outer Lindblad resonance (OLR) at r:!290" (6.3 kpc), 
Figs 21 and |5l The OLR corresponds to the scale of the 
optical and molecular gas disk. 

While introducing a low-mass nuclear component has 
a negligible effect on the rotation curve, it has a noticeable 
effect on the angular- velocity diagram. In the case with the 
single-disk mass distribution (solid line in the upper panel 
in Fig.^ there is no indication of an inner Lindblad reso- 
nance (ILR), which arises where il — k/2 = flp. However, 
the CO emission at the nucleus of M 83 is sph t into two 
distinct components i Petitoas fc WilsonI l| 19981) : Paper I), 
and it has been shown that the location (and existence) 
of a double nucleus of molecular gas in b arred galaxies, 
is coi ncident with the location of an ILR ijKennev et al.l 
In addition. lElmegreen et all l|l998^ observed a nu- 
clear ring in a (J — K) color image of M 83. The radius of 
this ring is 9" and the authors proposed that it coincides 
with the outer ILR (oILR). The disk plus central compo- 
nent model does allow for the existence of an ILR (lower 
panel in Fig. 2)). In the example above with an inner disk 
mass of 3 X 10^ Mq, the oILR would be located at a ra- 
dius of « 20". The possibility of an ILR is therefore not 
inconsistent with our data, but its exact location cannot 
be determined (nor can its existence with any certainty) 
due to the limited spatial resolution and grid spacing of 
our data. 

Figure 0] indicates the presence of an inner 
Ultraharmonic Resonance (iUHR), where f2— k/4 = fJp, at 
a radius of 60"±20" (1.3±0.4 kpc) (see also Fig.©. (The 
distance from the nucleus to the iUHR is only about 5" 
smaller if the low-mass nuclear component is not present.) 



In this radial range in the SW part of the bar a very 
conspicuous horn-shaped feature can be seen in both our 
C0( J=2-l) map (Fig.d) and in Ha images. In the NE sec- 
tion of the bar there is an abrupt decrease of the Ha emis- 
sion along the leading edge of the bar between the radii 
45-55", and in the C0(J=2-1) map the emission shows 
a sharp bend at this location. Furthermore, between radii 
40" and 60", the C0(J=2-1) intensity drops by a factor 
of about 4, and the corresponding drop in the C0(J=1- 
0) intensity is about 2-3. It is possible that these features 
are the results of the interaction between gas and stars on 
square-shaped 4:1 orbits and material flowi ng along the 
leading edge of the bar towards the nucleus. iRegan et alJ 
(1999) offer another explanation: their SPH- modeling sug- 
gests that some gas that initially follow the dust lane on 
the leading edge of the bar will over-shoot the nuclear re- 
gion and end up in the "spray region" on the opposite side 
of the bar in a region behind the dust lane. When this gas 
encounters gas flowing along the leading edge of the bar 
their trajectories will be more or less orthogonal to each 
other, and the interaction could have a strong impact on 
the gas flow in the dust lane, particularly if the gas feeding 
process is intermittent. 

We note that the spiral pattern of M83 appears to 
avoid the CR region. The spiral arms cross the CR region 
at (100", 150") and (-80", -130") in the coordinate system 
of Fig. |S1 The pitch angle increases abruptly when the 
arms approach the CR and decreases again on the other 
side of the CR region. Upon encountering the OLR, the 
pitch angle decreases and becomes negative. The spiral 
pattern shows here the characteristic dimpled shape of the 
2:1 orbits that run in the r etrograde sense in the rotating 
frame inside the OLR, e.g.. iKalnaisI l)l99l|) . 

The decrease of the angular velocity around 40" («1 
kpc) could be due to the effect of streaming motions in the 
bar. Gas on xi orbits, i.e. elongated orbits along the bar, 
will have a rotational velocity at the apex of the ellipse 
which is lower than that of gas on a circular orbit at the 
same radius. Hence, the ft will be lower. In the SW data, 
at galactocentric distances between 2 and 3 kpc, all mea- 
surements indicate deviations from circular motion. This 
can also be seen in the angular-velocity diagram: the NE 
sector shows a steady decrease of the angular velocity as a 
function of radius, while the SW sector has a region where 
the curve is more or less flat. This is the location of the 
western bar end, and we will discuss a possible explanation 
for this behavior in the next section. 



3.4. Isovelocity maps and residual velocities 

In Figs El and [3 the isovelocity contours obtained from 
the deconvolved CO(J=1-0) and C0(J=2-1) data sets, 
respectively, are plotted on grey-scale maps of the CO 
line emission in the respective transition. Both maps show 
the characteristic pattern of an inclined, rotating disk. 
Wiggles, indicating deviations from local circular motion, 
are present in the entire disk. Along the leading edge of 
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^Cdl. Distdnce; (kpc) 

Fig. 3. Circular velocities for the convolved C0( J=1~0) (upper panel) and MEM-deconvolved C0( J=2~l) (lower panel) 
data sets (within 40° of the major axis, binned into 6" (130 pc) wide annuli) shown as circles. Filled and unfilled circles 
represent data from the north-east and south-west sector, respectively, and the error bars are la. Circular velocities 
from Ha [stars] (de Vaucouleurs et al. 1983) and H i [triangles] (TA) are also shown. The solid lines are fits of a rotation 
curve generated by an exponential disk (upper panel) and a disk + central component (lower panel) model. 



the bar (see Fig.EJ), the shape of the iso- velocity contours 
points towards lower circular velocities than expected. 
This is caused by non-circular motions due to the oval 
gravitational potential in the bar and/or a shock asso- 
ciated with the leading edge of the bar. At the center 
of the galaxy the iso-velocity contours are tilted by 45° 
with respect to the minor axis, suggesting strong non- 
circular or out-of-plane motions (the former is expected 



in the presence of a strong bar). In CO w e do not detect 
any strong evidence of the warp seen in H i llR,og-stad et alJ 
'l974':'Huchtmcier fc Bohnenstengef 198lHTilanus & AllenI 
JJ93; Crosthwaite et al. 2002). This is not an unexpected 
result since a warp mainly affects regions outside the op- 
tical disk. 

Displayed in Fig. |H1 are the residual velocities, re- 
maining after subtraction of the axi-symmetric-model 
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Fig. 4. The angular velocity as a function of radius. The velocity data in the upper and lower panels are taken from 
the CO(J=1-0) convolved and the MEM-deconvolved C0(J=2~1) data sets, respectively. The solid line in the upper 
panel is the the best fit of the adopted rotation curve to the data (see Table EJ. In the case of the C0(J=2-1) data 
the rotation curve is composed of a disk and a central component (see Fig. Dashed and dotted lines represent: 
i7 — k/2, H. — k/4, si + k/4, and fl + k/2, where H. is the angular velocity and k the epicyclic frequency. The thick solid 
line in the bottom panel shows the adopted pattern speed (47 km s~^ kpc^^). 



velocity field, in the MEM-deconvolved CO(J=1-0) and 
C0(J=2-1) data sets, as well as the Hi data set from 
TA. In order to have comparable residual maps we sub- 
tracted the same axi-symmetric velocity field (the one ob- 
tained when fitting to the convolved C0( J=l-0) data set). 
Also shown in these figures is a line marking the crest of 
the gas mass surface density along the spiral pattern, as 



obtained from the CO(J=1-0) and C0(J=2-1) images. 
The residuals seem to follow a spiral-shaped pattern in all 
maps. In general, the residual velocities in the different 
maps show a good resemblance: both in the total range 
(from —20 to -1-20 kms~^) and in specific details, e.g., the 
large region with positive residuals 3f5 S of the center, 
and the bifurcation 2'.5 W of the center. The Hi map has 
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Fig. 5. The radii where our model predicts the locations of the inner ultra harmonic resonance (iUHR), corotation 
(CR), outer ultra harmonic resonance (oUHR), and outer Lindblad resonance (OLR) in M 83, on the velocity- integrated 
C0( J=l-0) intensity in the MEM-deconvolved data set. 

a slightly better resolution (12'.'4), but our maps show in 
addition the residuals in the central region where the H i 
emission is absent. Some regions in the residual velocity 
maps show very sharp gradients, such as the end of the 
arm bifurcation to the west. These areas can be identi- 
fied in the H I residual- velocity map as well. All of these 
regions coincide with areas in the disk with low mass sur- 
face density in both molecular and atomic hydrogen. Our 
interpretation is simply that the relatively low gas mass 



surface density in these regions is due to some mechanism, 
such as e.g. "shock focusing" (Lord & Kcnnev 199l'), that 
regulates how much gas can enter the region, and in the 
process skews the velocity distribution of the gas. 

An interesting aspect of the residual-velocity pattern 
is that it may be used to find the pattern speed. The 
technique is called "geometr ic-phase technique" , and it is 
described in ICan ziaia (^2^. According to this method, 
the residual- velocity pattern goes from being a one-armed 
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Fig. 6. Isovelocity curves (black lines) superposed on the velocity-integrated CO(J=1-0) intensity in the MEM- 
deconvolved data set (contour increment 10 kms~^). 



pattern to a three-armed pattern at the corotation radius 
in a two-armed spiral galaxy, where the pattern is driven 
by a spiral density wave (i.e., the pattern bifurcates at 
CR). On the western side there is such a bifurcation, but 
on the eastern side the situation is less clear. The galac- 
tocentric distance to the bifurcation on the western side 
is 160-180" (3.5-3.9 kpc), i.e., at 1.6-1.8 i?bar (i?bar is 
estimated as the radius where the isophotes in a K-band 
image become significantly influenced by the spiral arms) . 
This is inconsistent with the results from early analytical 



calculations ijContonoulos fc Grosbo]lll989l) . which show 
that corotation occurs just outside the end of the bar, and 
also Jarji;er than results from numerical SPH simulations 
l|Athanassoulalll99^ . which suggest a distance to the CR 
of 1.2±0.2 i?bar. Using Fig. ^we deduce that the angular 
pattern speed of M83 is about 40-55 kms~^ kpc""'^. This is 
in excellent agreement with th e value 50±9 kms~^kpc~^ 
derived by IZimmer & Rand j 12003^ by applying the 
Tremaine- Weinberg method ijTremaine fc WeinberE.198^ 
to our CO(J=1-0) data set. 
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Fig. 7. Isovelocity curves in the inner part of M83 superposed on the velocity-integrated C0(J=2-1) intensity in the 
MEM-deconvolved data set (contour increment 10 kms^^). 



3.5. Position-velocity diagrams 

The position-velocity (PV) diagrams are obtained along 
the position angles 45° (major axis, left row in Fig. EJ, 
and -45° (minor axis, right row in Fig. across the IR 
center of the galaxy. We used the deconvolved CO(J=1-0) 
and C0( J=2-l) data sets, but the behavior is identical in 
the convolved data sets, however, with a lower angular res- 
olution. We further assume in our interpretation that all 
motions take place in the plane of the galaxy (inclination 
24°), and that the NW side is the near side. 



3.5.1. Major axis 

Along the major axis of the galaxy the CO(J=1-0) PV- 
diagram crosses spiral arms twice, and there are velocity 
gradients over both arms. When compared to the expec- 
tations from density-wave-d riven streaming motions in a 
disk with trailing arms, e.g., (IVisseJlflSOaHCanzian 'lQQ?^. 
the effect is most pronounced in the SW arm. Outside CR 
and outside the arm, the gas is slowed down in the az- 
imuthal direction by the more rapidly moving spiral po- 
tential, while on the inside it is accelerated. In the PV- 
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Fig. 8. Residual-velocity maps after subtraction of the axi-symmetric component (described in the text) from the 

observed velocity fields in the deconvolved CO(J=1-0), deconvolved C0(J=2-1), and Hi data sets. The resolutions 
are 23", 13", and 12", respectively. The grey-scale and contour interval is —15 kms~^ (white) to +15 kms~^ (dark). 
The contour increment is 5 kms~^ and the thick contour represents kms~^. In the top image the position of the 
nucleus is shown as an encircled cross, and the positions of the bar ends are marked with circles. In all images the 
spiral arms, as seen in CO radio line data, are drawn with a dashed line. 
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diagram this appears as an increase in velocity, above that 
of the rotation curve, before the arm and a decrease in ve- 
locity, below that of the rotation curve, beyond the arm 
(in the direction of increasing galactocentric radius; note 
that the exact location of the spiral potential minimum is 
somewhat uncertain and not necessarily coincident with 
the CO arms which are indicated in the figure). The total 
azimuthal velocity shift is about 15 km s~^, or about 35 
km s~^ in the plane of the galaxy. This i s comparable to 
the streaming motions o bserved in M 81 l)VisseHh9Rnlil : 
lAdler fc WestpfahllllOQfil), but smaller than the ones ob- 
served in M51 jR ^bcck ct al. 1985: Garcia-Burillo et alJ 
ll993bHKuno fc Nakai,.1997; .Aalto et al...l9991 The situ- 
ation is less clear-cut in the NE arm, but the behavior is 
not inconsistent with that of density-wave-driven stream- 
ing motions. One explaination for the lack of expected 
streaming signatures at the NE crossing might be due to 
its proximity to the CR radius. Here the gas is comov- 
ing with the gravitational potential and thus subject to 
strongly non-linear effects. 

In the central region the velocity changes by about 50- 
75 kms~^ over «30" when crossing from the SW to the 
NE side of the nucleus (the effect is most pronounced in 
the C0(J=2-1) map). If this motion lies in the plane of 
the galaxy, the rotation speed is 60-90 kms~^ at a dis- 
tance of 15" (330 pc) from the nucleus. This suggests a 
much more rapidly rising rotation curve in the central re- 
gion than the one obtained from the azimuthally-averaged 
data in Sect. 13.21 and the presence of a central mass com- 
ponent (as discussed in Sects 13.21 and I3.3|l . Assuming a 
rotational speed of 60 kms~^, the enclosed mass within 
330 pc is about 3 x 10^ Mq. The emission in this region 
is dominated by the two centr al components observed i n 
e.g. CO(J=3-2 and J=4-3) bv lPetitr.a,s fc WilsonI ill 998l) . 
and in C0(J=2-1) by us (Paper I). Note that the veloc- 
ities of these concentrations, about 480 and 555 kms""'^ 
in the C0(J=2-1) map, are not symmetrical with respect 
to the systemic velocity, 14ys=510 kms^^). This can be 
explained by the fact that the two central components lie 
at different radial distances from the kinematical center. 
The result is also consistent with the CO(J=1-0) PV- 
diagram of lHanda et a l.' ('l990j), and th e hig her- J CO line 
PV-di agrams of llsrael fc Baaa ()200l|) and iDumke et alJ 
ll200l|) . 

In the lower left panel of Fig. we also show 
the rotation curve obtained using the envelope-tracing 
metho d to estima,te vel ocities in a major-axis PV dia- 
gram llSofue et a,l.lll997|> . The first step is to locate the 
terminal velocity (ft), i.e., the velocity at which the in- 
ten sity is 20% of the pea k intensity, for each radius (see 
e.g. lSofue fc RubinI l)200l|) for details of the method). The 
circular velocity {vc) is then calculated using 



= {vt - Wsys)/ sini - (0-2^, + CTi^sm)^^^ 



(2) 



where ctism is the velocity dispersion of the ISM, and CTobs 
is the velocity resolution of the observations. We used the 
radial distribution of the total velocity dispersion shown 



in Fig. ^1 as an estimate of (Cobs + ^ism)^^^- "^^^ 
cular velocities obtained in this way are systematically 
higher than the ones obtained with the intensity-weighted- 
velocity method. The shape of the rotation curve is con- 
sistent with a massive central comp onent. An even ste eper 
rotation curve has been derived bv lSofue et al ](|l99i) us- 
ing this method on H i and CO data. The main reason for 
the difference lies in the higher spatial resolution of the 
latter observations. Finally, it should be noted that this 
rotation curve is obtained along the bar, which happens 
to lie along the major axis. In this region there may also 
be non-negligable motions perpendicular to the plane. 



3.5.2. Minor axis 

Along the minor axis the CO(J=1-0) PV-diagram crosses 
the arms four times; denoted, starting from lower left, 
SEo (South-East outer), SEi (South-East inner), NWi, 
and NWo. 

In the SEi and SEo crossings the behavior is once again 
as expected from density- wave-driven streaming motions. 
The gas is accelerated radially outwards inside the arms, 
and radially inwards outside the arms. The velocity shifts 
are about 15 km s^^, corresponding to about 35 km s^^ 
in the plane of the galaxy, i.e., consistent with the velocity 
shift estimated at the major axis NE crossing. In the NW 
part the situation is much more complicated. The NWo 
crossing may be aff ected by a warp of the plane as sug- 
gested by Hi data llTilanus fc Allenlll993l) . Between the 
NWo and NWi crossings there is gas with radial velocities 
approaching 75 km in the plane of the galaxy, both in- 
wards and outwards, and this makes the interpretation of 
the NWi crossing more difficult. The gas moving inwards 
comes from an area slightly east of the minor axis, in a 
region where the arm bends sharply. Presumably this has 
an effect on the kinematics. In general, the gas is much 
more perturbed in the NW part. 

3.6. Velocity dispersion 

We estimate the velocity dispersion in the spectra by mea- 
suring their equivalent width, i.e., the velocity- integrated 
intensity divided by the peak intensity. This is converted 
to the standard, intensity-weighted, velocity dispersion, 
"■cO: by dividing the result with \/27r, i.e., the conversion 
constant for a Gaussian profile. Prior to the calculation 
of the velocity dispersions the spectra were smoothed to 
5 kms~^ in order to minimize the impact of the noise. 
This also facilitates the comparison of data with different 
velocity resolutions. 

Figure ^| shows the velocity dispersions of the 
CO(J=1-0) and C0(J=2-1) lines, at a resolution of 49", 
as a function of galactocentric radius. In the central re- 
gion the lines are markedly broadened by the steeply ris- 
ing rotation curve and bar-induced non-circular motions. 
In the disk there is a steady decrease of the velocity dis- 
persion with increasing galactocentric radius, reaching a 
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Fig. 9. Position-velocity diagrams along the major (left column) and minor (right column) axes, across the IR center, 
in CO(J=1-0) (upper panels) and C0(J=2-1) (lower panels) emission. The contours are 0.035, 0.07, 0.14, 0.3, 0.6, 1.2, 
1.8, 2.4 and 2.8 K (first white contour is 0.6 K). In order to show the details at the center, the y-scale is expanded in 
the C0(J=2-1) panels. The rotation curves obtained in Sect. are drawn as solid lines in the major axis diagrams. 
Additionally, in the lower left panel, the rotation curve obtained with the envelope-tracing method is also drawn. The 
latter rotation curve can be identified by its systematically higher higher rotational velocities. 



value of about lOkms^^ at the edge of the optical disk. 
The dispersion of the C0(J=2-1) lines is slightly higher 
than that of the CO(J=1-0) lines, by about 0.5-lkms~"'^ 
in the disk. This could possibly indicate the presence of 
several cloud components with slightly different physical 
characteristics. 

Before taking the decrease of the velocity dispersion 
in the disk with increasing galactocentric distance at face 
value, we need to evaluate the influence of the beam smear- 
ing on the velocity profiles. In the convolved CO(J=1-0) 
data set (upper left map in Fig. Ill|l . the pattern of high 
velocity dispersion is orthogonal to the bar, but in the 
deconvolved data set (upper right map in Fig. Ill|l the ef- 
fect is not present. This suggests that this feature is only 
an effect of the beam coupling to the intensity distribu- 
tion in a region where the local velocity and intensity 
changes within the beam are large. We made two sim- 
ple models to show the impact of this effect. In the first 
model we assumed a flat intensity distribution in order to 



look at large-scale patterns, and in the second we took 
the intensity distribution from the peak intensity distri- 
bution in our deconvolved data set in order to show arti- 
facts that may appear on smaller scales. In both models 
we assumed a fixed velocity dispersion of 8 km s^^ (see 
Sect. 01 and the velocity field generated by the rotation 
curve of the CO(J=1-0) convolved data. We created two 
data cubes, convolved them with a HPBW of 45" and 
calculated the velocity dispersion in the spectra using our 
automated routine (see Paper I for details). These maps of 
the model velocity dispersions can be seen in the bottom 
row of Fig. ^2 In both cases the velocity dispersion distri- 
bution is box-shaped and oriented with the "major axis" 
orthogonal to the bar. At the center the velocity dispersion 
reaches 20 kms~^ in both cases. The model spectra in the 
central region reproduce the shape of the observed spectra 
nicely. However, they underestimate observed velocity dis- 
persion at the center (40 km s~^), which is not surprising 
given our simple assumptions of the intensity distribution, 
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Fig. 10. Velocity dispersion in the C0( J=l-0) and C0( J=2-l) data sets, at a common spatial resolution of 49", as 
a function of the galactocentric distance. 



rotation curve, and low velocity dispersion in the nuclear 
region. In the model with the varying peak intensity, the 
elongation orthogonal to the bar is more pronounced. The 
explanation is that spectra in regions with low levels of 
emission, such as regions close to the bar (along the mi- 
nor axis), are influenced by the velocity components in 
nearby regions with high levels of intensity. This effect 
is also present in various regions of the disk. The lack of 
emission between the GMAs (introduced in Sect.0J and in 
the interarm regions leads here to an artificial broadening 
of the spectra. This results in iso-velocity-dispersion con- 
tours with fingers pointing outwards between the GMAs. 

The fact that the artificial broadening of the spectra 
decreases with galactocentric radius indicates that beam 
smearing is responsible for at least some of the gradient 
in the the velocity dispersion in Fig. ^| Our models did 
not include streaming motions, which also will add to the 
observed velocity dispersion. 

4. Disk instability 

Along the spiral arms we find clumpy structures, which 
we interpret as Galactic Molecular Associations (GMAs). 
Figure ^1 shows the locations of these structures on a 
grey-scale image of the total gas (H2+Hi+He) mass sur- 
face density (Sgas). The locations of the GMAs have been 
estimated from the velocity-integrated CO(J=I-0) emis- 
sion in the MEM-deconvolved data set. In some ambiguous 
cases we have used the peak intensity in the same data 
set. The II2 mass surface density is estimated from the 
velocity-integrated CO(J=1-0) line intensity using a con- 
version factor (Xco) of 2.3 xlO^° (Kknis^^)^^ cm~^ (see 
Paper I). The Hi data comes from TA. It was noted by 
TA that due to the lack of short-spacing information in 
the VLA ob servations only 55% of the at o mic g as mass 
observed bv iHiichtmeier fc Bohnenstengell l(l98ll) with a 
single dish telescope was detected. In order to calculate 



the total gas mass surface densities, we have to compen- 
sate for the missing flux. The difference in atomic masses 
detected in these observations corresponds to an average 
atomic mass surface density of 2.6 Mq pc^^. Lacking in- 
formation on how this mass is distributed in the galaxy, 
we take the simplest approach and add this constant mass 
surface density to the TA data. This is a minor correction, 
since the typical mass surface density in the GMAs are of 
the order 30 Mq pc~^. The final data are convolved and 
regridded to match the sampling of the CO data. The to- 
tal gas mass also includes He, through a scaling of the 
total Hi and H2 gas mass with 1.36. 

The total gas mass of an individual GMA is typi- 
cally 1-3 X 10'' Mq, and the distance between them are 
of the order 1 kpc. Similar structur es have been obser ved 
in other galaxies, such as e.g. M51 llKuno et allliool . It 
is believed that structures like these are formed when the 
gaseous disk becomes unstable to axisymmetric perturba- 
tions, which happens w hen the mass surface density ex - 
ceeds a critical density (lToomrelll964llKuTio et 3,1.11199,^1^ . 



where a is a dimensionless factor of order unity, Ugas the 
one-dimensional velocity dispersion of the ISM gas, and k 
the epicyclic frequency given by 

where 51 is the angular v elocity VrJR, and -Rp^ the c enter 
of the epicyclic motion ijBinnev fc Tremainelll987|) . For 
a pure stellar or gaseous disk one expects Q!=1, but in 
a realistic gas/stellar disk the value should be < 1 due 
to the instability intr oduced by the interaction between 
the two components llJoe fc Solomonlll984l: IWang fc Silk! 
I1994I) . lKenniciit!i] |l989*) studied the radial dependence of 
massive star formation and gas mass surface density in a 
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Fig. 11. Velocity dispersion in the CO(J=1~0) data (upper panels), two models showing the impact of the effect of 
the beam coupling to the rotation curve (bottom left) and in addition the intensity variation (bottom right). The 
contours range from 5 to 35 kms^^ and the thick line corresponds to 15 kms^^. The increment is 5 kms^^ in the 
upper panels and 2.5 kms^^ in the bottom panels (in order to enhance the effects in the models). 



number of galaxies and found that a lies in the range 0.5- 
0.85. An improved study bv lMartin fc Kennicut^l200l|) 
gave a = 0.69 (with some considerable scatter), assuming 
a gas velocity dispersion of 6 km s^^. In this paper we use 
a = 0.65. 



In a gravitationally unstable disk the GMAs are sepa- 
rated by a distance defined by the wavelength of the fastest 



growing mode in the gaseous disk l)Elmegreenlll994|) . This 
wavelength can be expressed as 



A = 2.2 



'gas 



7 km 



Jgas 



20 Mq pc- 



kpc. 



(5) 



Figure ^| shows the separation between the GMAs as a 
function of the average total gas mass surface density in 
the region. Also shown is the best fit of Q to the data. 
The fit gives an average velocity dispersion of 7.8 (±0.9) 
km s~^, only slightly lower than the values 9.3 km s~^ 
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Fig. 12. Map of the total gas (molccular+atomic) mass 
surface density in M83 with the position of the GMAs 
marked with squares and circles for the arms originating 
east and west of the nucleus, respectively. The contour 
values are 10, 20, 30, 40, 50, 75, 100, 200, 300, 400, and 
500 Mq pc-2 



found for the H i gas in TA, and 8.0 km s^^ for the CO in- 
terferometer data in lRand et alJ l)l999l) . and in agreement 
with the lower envelope of the values in Fig. II II This result 
als o gives some cred ence to our Egas estimates. 

lElmegreenl l|l994j) gives an expression for the character- 
istic mass of the GMAs as a function of velocity dispersion 
and mass surface density. 



20 Mq pc- 



M0.(6) 



This mass agrees well with our GMA mass estimates in 
M83, 1-3 X 10^ Mq. 

Using ||2J), setting tTgas = 7.8 km s~^, and using the ro- 
tation curve derived in Sect. 13 we have calculated the 
critical mass as a function of position in M83. Figure 
shows a map of the ratio T = Egas/Scr. T is below unity 
between the arms, hence these regions are stable against 
gravitational collapse. In the arms, T is usually about 1- 
2 and the variation with galactocentric distance is small, 
as expected since star formation drives T towards 1. The 
bar ends and the nucleus show values significantly larger 
than 1. These regions also have massive star formation, 
but the ratio is most likely overestimated since the ve- 
locity dispersion is higher than the adopted 7.8 km s^^ 
in these regions (especially close to the nucleus). Also, 
the true rotation curve is expected to be steeper than the 
one derived from our CO data, which would give a higher 
K and therefore an increase of the critical density. The 
value of T correlates nicely with the location of the H ii- 
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Fig. 13. The GMA separation as a function of mass sur- 
face density. The line shows the best fit of (0) to the data 
points. 



regio ns, which are shown as circl es in Fig. ^] (H ii data 



from lRiimstav fc Kaufmai] 



5. Conclusions 



)• 



We have used CO radio line emission data to investigate 
the kinematics of the molecular gas over the entire optical 
disk of M83. We compare our data with that of other 
kinematical tracers, such as Ha and atomic hydrogen. 
We find that: 

i) The large scale velocity field is regular. Deviations 
from circular rotation in the disk is caused by stream- 
ing motions associated with a spiral density wave. In the 
center region the velocity field is strongly coupled to the 
barred potential. 

ii) The kinematical center of the CO velocity field co- 
incides with the IR nucleus of the galaxy. The dynamical 
mass of the galaxy within the Holmberg radius is found 
to be 6 X 10^" Mq. The total gas mass, including II2, Hi 
and He, makes up 13% of this mass. 

Hi) The rotation curve is symmetric and adequately 
fitted by an exponential mass distribution. The rotation 
curve derived from the CO emission agrees with those es- 
timated from Ha and Hi emission. 

iv) Subtracting the global, rotational velocity field we 
find that the magnitude and pattern of the residual ve- 
locities of the CO line emission agree with those of H i. 
We also find that the CO residual velocity pattern bifur- 
cates on the western side. This is expected for a pattern 
driven by a spiral density wave with a co-rotation radius 
at the galactocentric distance of the bifurcation. However, 
no strong bifurcation is found on the eastern side. 

v) Position-velocity maps along the major and minor 
axes of M83 show velocity gradients across the spiral arms 
in agreement with density wave theory. The velocity shifts 
across the spiral arms are 35 km s""'^ (corrected for incli- 
nation) . 
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Fig. 14. A map of the ratio T = Sgas/Scr- Regions with T > 1 are unstable to axisymmetric perturbations and should 
collapse du to self- gravity of the gas. The contours represent T = 0.5, 1, 1.5, 2, 3, 4, and 5. The first white contour is 

T = 2. Also shown as circles arc the location of the Hii regions in M83. The size of a circle reflects the luminosity of 
the Hii region, ranging from non-giant Hil regions (small circles) to giant Hii regions (large circles). 



vi) The velocity dispersion in the disk appears to show 
a monotonic decline with galactocentric radii. Part, or 
all, of the decline can be attributed to a convolution of 
the velocity field and the finite telescope beam. The ve- 
locity dispersion of the C0(J=2-1) line emission is, how- 
ever, systematically larger than that of the CO(J=1-0) 
line emission by 0.5-1 km/s. 



vii) Several Giant Molecular Associations (GMAs) are 
identified along the spiral arms of M83. The separation of 
these GMAs depends on the velocity dispersion and the 
total gas surface density. Our derived CO velocity disper- 
sions and estimated gas siirface densities, together with 
the observed GMA separations, give a good fit to the ex- 
pected wavelength of the fastest growing mode of a grav- 
itationally unstable disk. 
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via) The total gas surface mass density, derived from 
CO and Hi measurements, is compared to the critical 
mass surface density as expected from a Toomrc disk. We 
find that the observed gas surface mass density exceeds 
the Toomre value in the spiral arm regions, but not in the 
interarm regions. We also find that the number and sizes 
of H II regions correlate well with the value of the ratio 

^gas/^cr* 
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